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ABSTRACT 

We have obtained 11.7 /irn and 17.9 /im images at the Keck I telescope of the 
circumstcllar dust emission from L 2 Pup, one of the nearest (D = 61 pc) mass¬ 
losing, pulsating, red giants that has a substantial infrared excess. We propose 
that the star is losing mass at a rate of ~3 x 10“' M 0 yr -1 . Given its relatively 
low luminosity (~ 1500 L 0 ), relatively high effective temperature (near 3400 K), 
relatively short period (~ 140 days), and the inferred gas outflow speed of 3.5 km 
s -1 , standard models for dust-driven mass loss do not apply. Instead, the wind 
may be driven by the stellar pulsations with radiation pressure on dust being 
relatively unimportant, as described in some recent calculations. L 2 Pup may 
serve as the prototype of this phase of stellar evolution where it could lose ~15% 
of its initial main sequence mass. 

Subject headings: circumstcllar matter -stars: mass loss 


1. INTRODUCTION 

Red giants with luminosities near 10 4 L 0 typically display infrared excesses that are 
the result of dust formed in a stellar wind reprocessing optical light emitted by the central 
star (Habing 1996). The mass loss rates can exceed ICC 5 M 0 yr~ 4 and the outflow speed is 
typically near 15 km s _1 . These winds are important in stellar evolution because the star 
may lose over half of its initial main sequence mass during this phase. 

A standard model to explain this mass loss is that as the stars pulsate, shock waves form 
in the atmosphere which drive matter above the nominal photosphere. In the post-shock 
region, the matter cools and grains form. After the dust is created, radiation pressure then 
expels the matter into the interstellar medium (see, for example, Laniers & Cassinelli 1999, 
Willson 2000). Although the standard model successfully explains the properties of many 
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red giants winds, there are some systems where it does not apply (see Jura & Kahane 1999). 
This “standard model” does not consider the effects of binaries, rotation or magnetic fields. 

In a recent set of calculations of models for mass loss from pulsating red giants, Winters 
et al. (2000) found a class of solutions, that are distinct from the standard model and which 
the authors denote as “case B”, where the mass loss rate is less than 3 x 10 - ' M 0 yr^ 1 at 
an outflow speed less than 5 km s -1 . In these models, mechanical energy from the pulsations 
drives the mass loss; radiation pressure on grains is a relatively minor effect. These “case B” 
models apply to stars of relatively low luminosity, typically where L*/M* < 3500 L 0 /M 0 , the 
pulsational periods are less than 300 days and the effective temperatures of the mass-losing 
stars are greater than 2600 K. 

Winters et al. (2000) were not able to identify any stars which clearly exhibit “B-class” 
mass loss. Here, we report detailed mid-infrared observations obtained at the 10m Keck 
telescope of L 2 Pup, a relatively nearby pulsating, mass-losing red giant. We propose that 
our data for the circumstellar shell around L 2 Pup can be best interpreted as a “B-class” 
wind. 


2. L 2 Pup 

L 2 Pup (= HD 56096 = HR 2748) has a distance measured by the Hipparcos satellite 
of 61 ± 5 pc, thus making it, along with R Dor, one of the pulsating, mass-losing red giants 
that is nearest to the Sun. In the Yale Bright Star Catalog it is classified as M5IIIe and 
assigned a visual magnitude of 5.10. It is a semiregular variable with a period of 140.6 days 
according to the General Catalog of Variable Stars, but with a 137 day period according to 
Whitelock et al. (2000). 

With a time-averaged measure that rri K = - 2.24 mag and thus m bo i = 0.73 mag (White- 
lock et al. 2000), and if M Bo i(Sun) = +4.72 mag, then the luminosity of L 2 Pup is 1500 L 0 . 
The time-averaged value of (J-K) is 1.31 mag (Whitelock et al. 2000), and therefore, from 
the models presented by Bessell et al. (1989) and Houdashelt et al. (2000), the time-averaged 
effective temperature of L 2 Pup is about 3400 K. With these values of the luminosity and 
temperature, the time-averaged radius of L 2 Pup is 7.8 x 10 12 cm. 

Since L 2 Pup has a velocity relative to the Galactic Plane of 90 km s 1 (Feast & 
Whitelock 2000), it probably belongs to the thick disk population, and thus the mass of the 
main sequence progenitor of the star was about 1 M 0 (Jura 1994). L 2 Pup’s location on 
the H-R diagram can be reproduced by interpolating from the Lattanzio’s (1991) models for 
the Asymptotic Red Giant Branch evolution of a star with a metallicity that is 0.5 solar. 
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In these calculations, the star has a luminosity between 1000 L 0 and 1500 L 0 for about 5 
x 10 5 yr. Given its luminosity and effective temperature, it is also possible that L 2 Pup 
is a first ascent red giant (see the tracks by Girardi et al. 2000); the luminosity and color, 
do not uniquely indicate the star’s evolutionary phase. The detection of Tc would indicate 
that the star almost certainly lies on the AGB; Lebzelter & Hron (1999) report a “possible” 
detection of Tc in the atmosphere of L 2 Pup. 

L 2 Pup is a visual binary with a companion which in the year 1913 was at 62" separation 
(Proust, Ochsenbein & Pettersen 1981). However, in the Hipparcos data, this companion 
now lies at an angular separation of 101" (Dommanget & Nys 2000). This much larger 
separation in the Hipparcos data implies that the two stars are not a physical pair. 

Evidence that L 2 Pup is losing mass comes from its infrared excess, the detection of 
circumstellar molecules, and its large, time-varying optical polarization. The value of Fj,(12 
pm)/F u 2.2 (/mi) is typically near 0.40 to 0.50, which is much larger than 0.07, the charac¬ 
teristic photospheric value for this ratio of M giants (see Jura, Webb & Kahane 2001). Since 
the circumstellar shell dominates the emission from the system for A > 5 /mi, the dust shell 
re-cmits approximately 10% of the total stellar luminosity. 

Detected circumstellar molecules include H 2 0, OH, SiO and CO. The line profiles of 
these molecules differ from each other because the H 2 0, OH, and SiO lines are produced as 
masers while the CO is thermally excited. The H 2 0 profile is variable with a velocity peak 
at Vlsr — 38 km s -1 and a second component that shows emission at velocities as low as 
Vlsr — 31 km s' 1 (Lepine, Paes de Barros & Gammon 1976, Gomez Balboa & Lepine 1986). 
The OH maser emission has been detected at Vlsr — 33-34 km s -1 (Dickinson et al. 1986). 
The SiO maser emission is variable and has been detected at Vlsr — 43 km s -1 (Balister 
et al. 1977) and as a triple-peaked source with velocities between Vlsr — 32 km s^ 1 and 
42 km s^ 1 (Haikala 1990). The CO emission peaks for the different rotational transitions 
have been reported between Vlsr — 32.9 km s -1 and 33.7 km s -1 (Kerschbaum et al. 1996, 
Kerschbaum & Olofsson 1999). The outflow velocity derived from the full extent of the CO 
emission has been measured to be between 2.0 km s -1 and 3.3 km s -1 . 

The SiO maser emission typically arise within 1-2 stellar radii of the photosphere (see 
Doeleman, Lonsdale & Greenhill 1998) and thus measures the motions induced by the stellar 
pulsations rather than the gas outflow. The velocities of the true circumstellar molecules 
range from 31 km s” 1 to 38 km s^ 1 ; the midpoint of this range of circumstellar gas velocities 
is 34.5 km s -1 . The extreme values of the circumstellar gas emission are determined from the 
H 2 0 masers, but the CO lines are nearly as wide as this spread in the H 2 0 radial velocities. 
Feast & Whitelock report that optical measurements of the star’s heliocentric radial velocity 
yield 53.0 km s' 1 which corresponds to Vlsr = 34.2 km s -1 . Within the uncertainties, it 
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appears that the center of mass velocity of L 2 Pup is 34.5 km s _1 and that circumstellar 
molecules move with speeds ± 3.5 km s -1 with respect to this central velocity. Although the 
outflow is not spherical, we have no information on the spatial variations of the gas outflow 
speed and we therefore assign an outflow speed of 3.5 km s _1 to the wind. This wind speed 
is much lower than the typical value of >10 km s -1 found for most mass-losing red giant 
stars (Loup et al. 1993). To illustrate the point that L 2 Pup has an unusually low wind 
velocity, we show in Figure 1 a histogram of the outflow velocities measured from the line 
widths of the circumstellar CO emission from a survey of semiregular and irregular variables 
(Kerschbaum & Olofsson 1999). In this survey, L 2 Pup has the narrowest line, although 
there are some other stars such as EP Aqr with narrow spikes on top of a broader profile 
(see also Knapp, Young, Lee & Jorissen 1998). 

At times, the optical light from L 2 Pup is over 10% polarized (Magalhaes et al. 1986). 
Since L 2 Pup is only 61 pc from the Sun, this large polarization cannot be produced by 
interstellar grains but must be intrinsic to the star and its circumstellar envelope. The posi¬ 
tion angle varies with time and has ranged from 156° through 0° up to 60° (Magalhaes et al. 
1986). The large optical polarization implies that the circumstellar envelope is anisotropic. 


3. OBSERVATIONS 

Our data were obtained on 2001 Feb 05 (LIT) at the Keck I telescope using the Long 
Wavelength Spectrometer (LWS) which was built by a team led by B. Jones and is described 
on the Keck web page. The LWS is a 128 x 128 SiAs BIB array with a pixel scale at the 
Keck telescope of O'/OS and a total held of view of 1072 x 10"2. We used the “chop-nod” 
mode of observing, and 2 different Liters centered at 11.7 fin 1 and 17.9 fin 1 with widths of 
1.0 fim and 2.0 fix n, respectively. Following Chen & Jura (2001), we used Capella (= HR 
1708) for flux and point-spread-function calibrations. For Capella the FWHM of the image 
was 0747 and 0749 at 11.7 fin 1 and 17.9 fa n, respectively. 

The images of L 2 Pup in the 11.7 fan in the 17.9 /an Liters are presented in Figures 2 
and 3. The 11.7 fan image of Capella is shown in Figure 4. Because the calibrator was much 
more compact than L 2 Pup, deconvolution of the data is not required for oLsets greater 
than 076 from the star since at these oLsets, at most azimuths, the central point source 
contributes less than 20% of the observed intensity. In the mid-infrared, L 2 Pup is brighter 
than Capella, and we cannot be sure that some of the very faint extended emission that 
we detected is real rather than an artifact. We therefore only show emission at positions 
where the observed intensity is at least 2% of the peak intensity. At both wavelengths, we 
see an extended morphology oriented along an axis at position angle 135° with an additional 
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distinct blob at position angle of 225°. 

Integrated over the image, we measure F„(11.7 /jrn) = 2500 Jy; this is comparable to 
the non-color corrected value of F !/ (12 /irn) = 2400 Jy measured with IRAS and the results 
from DIRBE (see Figure 7 below) that F I/ (12 /mi) ranges between about 2000 Jy and 2200 

Jy- 

If we assume the same atmospheric extinction for L 2 Pup and Capclla, our data imply 
that F I/ (17.9 fin i) = 1100 Jy. However, because of its southern declination, there was an 
air mass of 2.32 during our observations of L 2 Pup. Therefore, the atmospheric extinction 
correction for the flux at 17.9 /irn of L 2 Pup was substantial and the absolute level of our 
measurements at 17.9 /irn is uncertain. The data from the DIRBE satellite shown below in 
Figure 7 shows that F*,(12 /im)/F i/ (25 /irn) remains constant to within 5% over a pulsational 
cycle. Interpolating the DIRBE data between 12 /irn and 25 /irn indicates an average flux 
at 17.9 /.irn of 1700 Jy instead of our measured value of 1100 Jy. We therefore apply an 
uncertain correction factor of 1.5 to our measured values of the intensities at 17.9 /irn. 


4. MODELS 

To model the infrared maps of the envelope, we follow a standard prescription where 
the dust grains are heated by light from the central star and then re-radiate in the infrared 
(Sopka et al. 1984). Such models require estimates for the size and composition of the 
grains. Here, we use simplified semi-analytic models to show how different assumptions and 
parameters propagate into our derived results. 


4.1. Grain Properties 

An important parameter is the grain size. Daniel (1982) has modeled the optical polar¬ 
ization around L 2 Pup by spherical grains with radius a = 0.26 /irn with an assumed simple 
bipolar configuration of the circumstellar envelope. In fact, as shown by our infrared imag¬ 
ing, there is a blob at position angle 225° in addition to the extended bipolar morphology 
extended along an axis at position angle 135°. The spatial structure of the circumstellar 
dust therefore is more complex than presumed in Daniel’s models. Everything else being 
equal, the more complex geometry will result in less integrated polarization. Therefore, in 
order to reproduce the result that the polarization in the integrated light approaches 10%, 
the grains must be more efficient at producing polarization and therefore smaller than those 
postulated by Daniel. Jura (1996) suggested that the typical spherical particle around a 
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mass-losing oxygen-rich star has a radius smaller than 0.15 /mi. Here, we assume that most 
of the particles around L 2 Pup are similar to those around other oxygen-rich stars and thus 
for most of the particles a < 0.1 fin i. Since T* = 3400 K, this estimate of the grain size is 
such that a < A/ (2n) for all wavelengths of interest both for emission in the infrared and for 
absorption of the light from the photosphere. In this case, the temperature of the grains is 
insensitive to their size. 


A model of the grain emissivity is necessary to explain our infrared data. Unfortunately, 
the spectrometer on the IRAS satellite did not acquire a good infrared spectrum of this 
object (Sloan & Price 1998, Volk & Cohen 1989), and ground-based infrared spectra have 
not been reported in the literature surveyed with SIMBAD. Therefore, the magnitude of 
silicate features in the infrared dust emission cannot be easily quantified. If we define Q a b s 
as the ratio of the absorption cross section to the geometric cross section (see, for example, 
Spitzer 1978), then the relationship between the opacity, (cm 2 g' 1 ) and Q a bs is that 
= (3 <3i/)/(4 p s a) where p s (g cm -3 ) is the density of the solid grains which we take to equal 
3.3 g cm -3 (Kim, Martin & Hendry 1994). A simplified version of the circumstellar silicates 
described by Ossenkopf, Henning & Mathis (1992) and David & Pegourie (1995) can be 
expressed as: Q a b s = K\ v +1 for A < 9 fim and Q a b s = K 2 i ,+1 for A > 9 fim with K { /K 2 
~ 0.2. This simplified model for Q v implies values of yy from which the grain temperature 
can be inferred, although this model does not include any spectral features such as those 
produced by silicates. 

Given yy, then the grain temperature profile can be found following the prescription 
that the particles are heated by starlight and re-radiate in the infrared. For L 2 Pup with T e ff 
= 3400 K, then typically the heating occurs near A ~ 1.5 fim, while most of the re-emission 
in the infrared occurs for A > 9 fim. Therefore, following Sopka et al. (1984), the grain 
temperature, T gr , as a function of distance, R, from the star is given as: 


T 

± gr 


>*0.2t-, (R*\0A 
4AV * [ R ] 


( 1 ) 


A validation of this approach is to compare our simple model with more sophisticated calcu¬ 
lations by other authors. For example, for a star with T e ff — 3400 K, at D = 20 /?*, equation 
(1) yields T gr = 563 K. With detailed calculations for small silicate grains, Lorenz-Martins 
& Pompeia (2000) compute that at R = 20 /?*, T gr = 567 K. The closeness of this agreement 
is fortuitous, but it does allow that our simple model might be a useful approximation. 


We estimate the intensity of the circumstellar emission with the assumptions that the 
mass loss rate, M gr , and the grain outflow velocity, V gr , are constant. If we observe along a 
ray where the measured distance is denoted by z, then the observed surface brightness, /„, 
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(Jy ster x ) is: 

! v = J XuP(z)B u (T(z))dz (2) 

For a spherically symmetric distribution, 


P(z) 


Mgr 

4 7T (6 2 + z 2 ) Vgr 


( 3 ) 


where b denotes the impact parameter of the sight-line with respect to the central illuminating 
star. The solution to equations (2) and (3) for M gr is: 


Mg r 


4 b Vgr Iu fcorr (jz ) 
Xu B u [Tgr(b )] 


( 4 ) 


The function f corr ( v) is a measure of the effects of a non-uniform temperature distribution of 
the dust grains and must be computed numerically. If the grains had a uniform temperature, 
then fcorriy ) — 1- For L 2 Pup, we estimate that at V from the star that T gr = 280 K. In 
this case, / corr (11.7 /mi) = 2.0 while / corr (17.9 /jrn) = 1.7. 


If Xu is known, then we can invert equation (4) to find from our measurements of l u , 
the value of T gr as a function of b. Assuming some silicate emission, we expect from Draine 
(1985) and Ossenkopf et al. (1992) that x„(11.7 y um)/x iy (17.9 /mi) = 1.6, instead of 1.5 as 
predicted if Xu oc v. With this ratio of the opacities and with the measurement that A (11.7 
/nn)/A(17.9 /mi) = 2 ± 0.7 at offset 1" from the star, we derive that at this location, T gr 
> 350 K. From equation (1), we predict that T gr = 280 K. Given the uncertainty in the 
calibration of the intensity at 17.9 /mi, and our lack of constraints on the grain composition, 
the simplified model and the data are consistent with each other. 

To estimate M gr , we need to infer the absolute value of Xu as well as its variation as 
a function of frequency. From Ossenkopf et al. (1992), we adopt x^(11.7 /mi) = 2000 cm 2 
g _1 . In the simplified model for the particle cmissivity described above, at the peak of the 
photospheric emission, Xu( 1-5 /mi) ~ 3000 cm 2 gW 


4.2. Mass Loss Rate 


To estimate the mass loss rate from equation (4), we must determine V gr . In the 
envelopes of red giants the grains stream supersonically through the gas, and the grain 
velocity is therefore larger than the gas outflow velocity, V ga . If the gas outflow rate is M ga , 
and if a < A/(27r) so that scattering by the grains is unimportant, then (Laniers & Cassinclli 
1999) 


V g r Fga 


- Qabsl-'* V ga \ 1/2 


M ga c 


( 5 ) 



Assuming the emission from the star with T e ff = 3400 K peaks near 1.5 /irn, and with Xu 
(1.5 /irn) = 3000 cm 2 g -1 , V ga = 3.5 km s -1 , M ga = 3 x 1O -7 M 0 yr -1 (see below) and a = 
10 -5 cm implying that Q a bs = 0.13, then V gr ~ 10.0 km s -1 . 

We measure at 1” from the star that A (11.7 /.irn) ranges from 140 to 280 Jy arcsec -2 . 
Adopting a “typical” value of 200 Jy arcsec -2 , and with V gr = 10.0 km s -1 , \v (H-7 /rm) = 
2000 cm 2 g -1 , and f corr ~ 2.0, then, from equation (4), M gr ~ 1.6 x 10 -9 M 0 yr -1 . 

We test this model by comparing its predictions with observations. As shown in Figure 
5, the model and the spatial variations of the observed intensity at both 11.7 fin l agree 
reasonably well with each other. However, as shown in Figure 6, the agreement between the 
model and the observations at 17.9 fim do not agree especially well. As discussed above, the 
data indicated a temperature at 1" of 400 K while the model predicted a temperature of 280 
K. Since the intensities at 17.9 /irn may be systematically in error because of the difficulty 
of calibrating our observations, we scale all the 17.9 /irn intensities upward by a factor of 
2 and plot these results in Figure 6 as well. In this case, the model agrees reasonably well 
with the observations. 

In addition to determining the dust loss rate, we would like to estimate the gas loss 
rate. The most common tracer of the amount of circumstellar gas is CO. However, the 
circumstellar envelope around L 2 Pup is unusual. Kerschbaum & Olofsson (1999) report 
that the integrated intensity in the CO (2-1) lines compared to that in the CO (1-0) line 
is 12, substantially larger than the value for this ratio measured for other semiregular stars 
and also larger than can be produced by any standard model for thermal excitation of this 
molecule (see Kahane & Jura 1994). Since we have only an uncertain measure of the amount 
of circumstellar gas, we adopt M ga /M gr « 200, consistent with the inferred evolutionary 
status of the star as a member of the thick disk. With a dust loss rate of 1.6 x 10 -9 M 0 
yr -1 , then the total mass loss rate is ~3 x 10 ' M 0 yr 1 . 

5. PULSATIONAL MODE 

L 2 Pup is a variable star, and the models of Winters et al. (2000) are consistent with 
our inferred mass loss rate. Here, we note that the energy carried away by the wind, 
(1/2) M ga I/ 2 , is about 1 x 10 3 ° erg s -1 . Wood & Karovska (2000) observed the Mg II 
ultraviolet emission from L 2 Pup during different phases of its pulsational cycle. This line is 
excited by shocks caused by the pulsations, and while variable, a typical flux at the Earth 
is 5 x 10 -13 erg s -1 . Therefore, shocks in the atmosphere radiate at least 2 x 10 29 erg s -1 , 
and thus there appears to be enough kinetic energy in the pulsations to drive the outflow, 



9 


lending support to the using the models of Winters et al. (2000) for explaining the wind. 

Red giants often exhibit asymmetrical mass loss (see, for example, Monnier, Tuthill & 
Danchi 2000, Uitenbroek, Dupree & Gililand 1998) whose origin is not well understood. It 
is possible that non-radial pulsations at least partly contribute to anisotropic mass loss, and 
this may be occuring around L 2 Pup. Here, we argue that the hypothesis of non-radial pul¬ 
sations is supported from observations of the time-variability of L 2 Pup at different infrared 
wavelengths. The basic idea is that the time variations in the flux from the circumstellar 
shell do not track very well the time variations in the flux emitted by the photosphere. One 
possible explanation for this difference is that the star undergoes non-radial pulsations. 

A comprehensive data base of infrared emission from bright stars was acquired with 
the DIRBE instrument on the COBE satellite where fluxes at 9 different infrared bands 
were measured, sometimes more often than daily, over a 300 day period. We have used the 
standard web-based tool to extract these data for L 2 Pup, and the results for 7 bands are 
shown in Figure 7. It can be seen that the source is substantially more variable at A < 3.5 
/mi, the light emitted by the photosphere, than at A > 4.9 /mi, the light emitted by the 
circumstellar dust. To illustrate the effect somewhat differently, we show in Figure 8, a plot 
of the deviation from the mean value of the DIRBE fluxes at three wavelengths, 1.25 /mi, 
2.2 /mi and 12 /mi. The fluxes at 1.25 /mi and 2.2 /mi track each other, while the flux at 12 
/mi displays a distinctly different time variation. 

The data in Figure 7 can be compared to a simple model for a spherically symmetric 
envelope where the specific luminosity of the dust, L u , from the optically thin envelope is: 

r°° M 

XM* B u [T gr ]){-f)dR ( 6 ) 

J 0 v gr 

If the temperature variation of the circumstellar dust is given by equation (1), then we can 
show from equation (6) with a simple substitution of variables that: 

L v oc Tf- 5 R* oc (L* T*) 1 / 2 (7) 


Consider now the observations of L 2 Pup shown in Figures 7 and 8. From the model 
atmospheres computed by Houdashelt et al. (2000), near T e ff — 3400 K, we expect that 


AT 

~T~ 


0.4 


Ar • 

L - Xl n.i.r 


( 8 ) 


where r nir denotes the flux ratio, Fj,(1.25 /mi)/F Jy (2.2 /mi). For example, from day 120 to 
day 190, the flux ratio, r n i r increases by a factor of about 1.08 which can be explained by an 
increase in the temperature by 1.03. Although the variation in radius during this phase is not 
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certain, in view of the total increase in the observed values of F„(1.25 //rn) and F^(2.2 /irn) 
as well as in r nir , it is probable that during this time interval, the radius slightly increased as 
well. Consequently, using equation (7), we expect from the inferred increase in T* that F 1/ (12 
/irn) should have risen by at least 8%. In fact, as can be seen in Figure 8, F^(12 /irn) at day 
190 is within 1% of its value at day 120. Furthermore, the time variation of the flux at 12 /irn 
is different from that at 2.2 /irn. Therefore, a model with radial pulsations is not supported 
by these DIRBE observations. One way to understand these data is to presume that the 
star undergoes nonradial pulsations. In this case, F^(12 /irn) measures the luminosity of the 
entire star while F^(2,2 /irn) only measures emission from the hemisphere of the star that 
faces the Earth. A test of this model is that the resolved dust blobs shown in Figures 2 and 
3 should display unsynchronized variations. 

Additional evidence for nonradial pulsations comes from the marked time variation of 
the position angle of the net polarization of the star (Magalhaes et al. 1986). These data 
can be explained if the circumstellar blobs are differentially illuminated by a non-spherical 
time-varying photosphere as would occur during nonradial pulsations. 

The DIRBE data can be used to study time-variations from other pulsating red giants 
besides L 2 Pup. We show, for example, in Figures 9 and 10, the DIRBE fluxes and their 
deviations from their mean values, for R Cas, a mass-losing infrared-bright Mira red giant 
with a period of 430 days which lies at a distance of about 110 pc from the Sun. This star 
exhibits time variations at both near infrared and mid-infrared wavelengths which follow 
each other as would be expected from radial pulsations. 


6. DISCUSSION 


With M ga « 3 x 10~ 7 M 0 yr -1 , V gas ~ 3.5 km s' 1 and L 1500 L 0 , the circumstellar 
envelope around L 2 Pup is well described by the models computed by Winters et al. (2000) 
where radiation pressure on the grains is relatively unimportant in the dynamics. To assess 
this hypothesis, consider a which denotes the ratio of the outward radiation pressure on the 
circumstellar envelope (gas plus dust) compared to the inward force of gravity. We may 
write (Winters et al. 2000) 


a 


MgrVgas / Ljf \ 

MgaVgr 47 TCGM* 


( 9 ) 


With the parameters given above, then for the envelope around L 2 Pup, a = 0.6 and radiation 
pressure does not dominate the circumstellar dynamics. 


Besides L 2 Pup, there are a few other evolved red giant stars where the measured 
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circumstellar gas exhibits velocity widths smaller than 5 km s” 1 . Jura & Kahane (1999) 
have proposed that some of these stars, such as the Red Rectangle and AC Her which are 
known close binaries (van Winckel et al. 1998), possess orbiting molecular reservoirs. They 
hypothesize that the orbiting material may be produced by interactions of the secondary 
star with the envelope ejected by the mass-losing star. There is no evidence for a companion 
around L 2 Pup. However, at 61 AU from the star, the region where our infrared data 
indicates the presence of a large amount of circumstellar dust, the escape velocity is 5.4 km 
s” 1 which is larger than the inferred gas outflow speed. If L 2 Pup has an undetected low-mass 
companion, the interaction of the wind with this star might produce some gravitationally 
bound matter in a circumbinary disk. 

An red giant star may be spun-up by engulfing a companion, and Soker (2000) has 
proposed that such stars with relatively high rotation rates may sometimes eject mass in a 
very flattened slow-moving outflow. However, there is no reported evidence that L 2 Pup is 
rotating particularly rapidly. Also, although this picture of a flattened outflow may partially 
explain the observation around L 2 Pup of an apparent axis of symmetry at position angle 
135°, the presence of a distinct blob at position angle 225° is not naturally explained with 
this model. 

The current mass-loss phase of L 2 Pup may be important in its evolution. If L 2 Pup 
persists in losing mass at ~ 3 x 10”' M 0 yr” 1 for 5 x 10 5 yr, consistent with the models 
for stellar evolution computed by Lattanzio (1991) of the AGB or by Girardi et al. (2000) 
for the tip of the Red Giant Branch, then it will lose ~0.15 M 0 or ~15% of its initial main 
sequence mass. This estimate of the lifetime of L 2 Pup in its current phase is completely 
theoretical since it takes the dust only ~ 100 yr to reach the projected separation of 2" from 
the star seen in our infrared images. 

L 2 Pup may serve as a test of the proposed mass loss rate from red giants. This is 
important because it is not possible to compute with confidence the mass loss rate from 
our current understanding of stellar evolution (see, for example, Catelan 2000, Schroder & 
Sedlmayr 2001). Although recognizing that there are many uncertainties, Catelan (2000) has 
proposed that the mass loss rate from red giant stars can be parameterized by the expression 

w) Meyr ~' (10) 

where g is the acceleration of the star (cgs), and L* and R* are in solar units. In the case 
of L 2 Pup, this formula predicts a mass loss rate of 1.1 x 10” 8 M 0 yr” 1 , a factor of 30 less 
than we infer. It could be that the formula of Catelan (2000) fails to describe adequately the 
mass loss from L 2 Pup because it does not distinguish between pulsating and non-pulsating 
stars. 


M ga = 8.5x10 


-10 



12 


Some red giants in 47 Tuc with roughly the same luminosity as L 2 Pup do not exhibit 
any 12 /mi excess (Ramdani & Jorissen 2001). A similar result is found for red giants in 
the neighborhood of the Sun (Jura et al. 2001). If these stars are losing mass, it is without 
the production of large amounts of dust. In the models of Winters et al. (2000), the dust 
is not necessary to drive the mass loss, and thus there might be red giants with large mass 
loss rates and little dust formation. Catelan’s formula possibly applies to non-pulsating red 
giants with luminosities near 1000 L 0 . 


7. CONCLUSIONS 

We have obtained mid-IR images of the dust around L 2 Pup. We find the following: 

• We estimate a total mass loss rate of ~3 x 10~ 7 M 0 yr _1 . During the current phase 
of its evolution, L 2 Pup may lose ~15% of its initial main sequence mass. 

• Non-radial pulsations may at least in part account for the observed asymmetric mass 
loss, time variations of the infrared fluxes, and time variations of the position angle of 
the optical polarization. 

• L 2 Pup might serve as the prototype of an outflow driven by stellar pulsations with 
radiation pressure on dust being relatively unimportant, as in models described by 
Winters et al. (2000). 

This work has been partly supported by NASA. 
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FIGURE CAPTIONS 


Fig. 1. Histogram of outflow velocity measured from the profile of the circumstellar CO 
emission for a sample of semiregular and irregular variables (Kerschbaum & Olofsson 1999). 
In the case of several measurements of the line profile, we adopt the average value. In the 
case of profiles where both a narrow and broad component were identified, we consider only 
the broad component. The location of L 2 Pup in this histogram is indicated. 

Fig. 2. The 11.7 /mi image of L 2 Pup. North is up and East is to the left. The con¬ 
tour levels (Jy arcsec -2 ) are shown in the color bar. 

Fig. 3. The 17.9 /mi image of L 2 Pup with the same conventions as Fig. 2. The very 
faint blob located at 2" east of the star is an artifact of the LWS filter at this wavelength. 

Fig. 4. The 11.7 /jrn image of Capclla with the same conventions as Fig. 2. 

Fig. 5. The measurements of 7^(11.7 /irn) at offsets > 0764 from the star at position 
angle 0° (circles) and 120° (crosses). The solid line shows the predictions for 7^(11.7 /mi) for 
the spherical model described in the text with M gr = 1.6 x 10“ 9 M 0 yr _1 , V gr = 10.0 km 
s^ 1 and T gr = 280 K at 61 AU from the star. 

Fig. 6. The same as Figure 5 except that we compare data at 17.9 /mi with the model 
shown as the solid line. As described in the text, the calibration of the data is uncertain; 
the red points show the same data scaled upwards by a factor of 2. 

Fig 7. The fluxes acquired by the DIRBE instrument on COBE for L 2 Pup. The labels 
refer to fluxes at different bands and are reported at different days during the COBE mis¬ 
sion relative to January 1, 1990. 

Fig. 8. The flux ratios measured by DIRBE for L 2 Pup relative to the mean value at 
different times for 1.25 /mi, 2.2 /mi and 12 /mi. The curves at 1.25 /mi and 2.2 /mi follow 
each other and track the emission from the photosphere while the 12 /mi emission is mainly 
produced by circumstellar dust. 

Fig. 9. The fluxes acquired by the DIRBE instrument on COBE for R Cas with the same 
conventions as in Figure 7. The time variations of the fluxes at 12 /mi and 25 /mi follow 
those at shorter wavelengths. 
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Fig. 10. The flux ratios measured by DIRBE for R Cas relative to the mean value at 
different times for 1.25 /nrl, 2.2 /mi and 12 /mi. 
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